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Abstract 

The cosmic-ray propagation code GALPROP has been generalized to include fragmentation networks of ar- 
bitrary complexity. The code can now provide an alternative to leaky-box calculations for full isotopic abun- 
dance calculations and has the advantage of including the spatial dimension which is essential for radioactive 

-^ . nuclei. Prehminary predictions for sub-Fe/Fe, ^°Be/ ^Be, ^^Al/ ^^Al, '^^Cl/Cl, and ^'^Mn/Mn are presented 

Qs ' in anticipation of new experimental isotopic data. 

o\' 

1 Introduction: 

The GALPROP program (Strong and Moskalenko 1998, Moskalenko and Strong 1998, Strong et al. 1999) 
performs cosmic-ray propagation calculations for nuclei, electrons and positrons and computes 7-ray and syn- 
chrotron emission in the same framework. The 3D spatial approach with a realistic distribution of interstellar 
gas distinguishes it from leaky-box calculations. Originally we considered that the complexity of the spatial 
model would preclude using cosmic -ray reaction networks with more than a few species, and our results so far 
^ ■ have hence been obtained using a simplifying 'weighted cross sections' method. However we have now been 

^ . able to generalize the scheme to include reaction networks of arbitrary complexity. The scheme can hence 

^ \ potentially compete with leaky-box calculations for computation of isotopic abundances, while retaining the 

\C> ' spatial component essential for radioactive nuclei, electrons and gamma rays. Even for stable nuclei it has the 

^ ■ advantage of a physically-based propagation scheme with a spatial distribution of sources rather than an ad 

Qv^ ! hoc path length distribution. It also facilitates tests of reacceleration. Here we give a progress report on the 

new extension of the method and present some illustrative results. The technique should be of interest for the 
pi i' interpretation of new measurements such as those becoming available from ACE and ISOMAX. 

g : 2 Method: 

^ . GALPROP solves the Galactic CR propagation equation numerically on a grid in 3D with cylindrical 

symmetry. In the original simplified approach for nuclei we propagated first the primaries, then the source 

function for each secondary was computed and then the secondaries were propagated. Multiple progenitors 

^ . and tertiary etc. reactions were handled by weighting the cross-sections according to the observed abundances, 

which was clearly a compromise solution. 

The new extended scheme instead handles the reaction network explicitly, as follows: 

1. Propagate all primary species from an assumed set of source abundances; 

2. Compute the resulting spallation source function for all species; 

3. Propagate all species using the sum of primary and spallation sources; 

4. Iterate steps 2 and 3 until converged. 
After the second iteration the result is already accurate for the pure secondary component, after the third 

iteration it is accurate for tertiaries, and so on. Hence in practice only 3 or 4 such iterations are necessary for a 
complete solution of the network. The method is more time-consuming than the simple approach since many 
more species are included and because of the several iterations required. 

3 Results: 

We have applied the method to a network of 87 nuclei from protons to Ni, including explicitly all sta- 
ble species and radioactive species with half-life more than 10^ years. The network of channels for decayed 
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cross-sections is based on Nuclear Data Sheets. Isotopic cross sections are based on measured values where 
available, up to and including Webber et al. (1998). For cases where measurements at enough energies exist we 
interpolate a smooth function. Otherwise we use the Webber et al. (1990) cross-section code, renormalizing 
to measurements where they exist. The interstellar He abundance is taken as 0.11 by number. Source abun- 
dances are from DuVernois and Thayer (1996), with solar isotopic ratios within a given element. Propagation 
parameters for this model with reacceleration are: halo size z^ = 4: kpc, va = 20 km s~^, diffusion coeffient 
D = DQf3{p/po)^, where Dq = e.TSxlO^*^ cm^ s~^, po = 3 GeV. The diffusion coefficient, adjusted to fit 
B/C, differs slighly from that used in the original work for the same Zh (Strong et al. 1998: Dq = 6x 10^^ cm^ 
s~^), reflecting the more detailed treatment and updated cross sections. 

Fig- n shows computed fluxes of all the included isotopes at 2 GeV/nucleon. This result is illustrative of 
the method but not to be taken as predictions for evaluation purposes. 

Fig. El shows B/C and sub-Fe/Fe. This model with reacceleration reproduces the sub-Fe data reasonably 
well considering that the model was adjusted only to fit B/C. The deviation from the observed sub-Fe/Fe 
shape is however noticeable, as discussed by Webber (1997). Figs. |3l and |4] show the radioactive isotope 
ratios ^°Be/ ^Be, 26a1/ ^^aI, 36ci/ CI, and ^^Mn/ Mn (using an estimated ^''Mn halflife of 6.3 x 10^ years, 
Wuosmaa et al. 1998). Modulation is for nominal values of the modulation parameter but this has little effect 
on the comparison due to the small energy dependence of these 4 ratios in the 100-1000 MeV/nucleon range. 
In Strong and Moskalenko (1998) we obtained a range 4 kpc< z^ < 12 kpc for the halo height; in the present 
improved model the ^'^Be/ ^Be, ^^Al/ ^^Al predictions are consistent with this. A complete grid of models 
will be generated in future to obtain new limits. Our ^^Cl/ CI and ^^Mn/ Mn predictions are preliminary and 
do not yet provide significant constraints. 

4 Conclusions: 

The extension of GALPROP to complete networks of nuclei has been accomplished, but further work is 
needed especially on the cross-sections before the method can be used for detailed evaluation. As usual the 
software will be made available on the WWW at http://www.g amma.mpe-garching.mpg.de/'^aws/aws.html| 
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Figure 1: Model isotopic fluxes at 2 GeV/nucleon as function of Z; individual A values are marked. NB: This 
plot is illustrative of the method but not to be taken as predictions for evaluation purposes. 
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Figure 2: Left: B/C interstellar and modulated to 500 MV for diffusive reacceleration model with Zfi = 4 kpc. 
Data compilation: Webber et al. (1996). Right: The same for sub-Fe/Fe (Z = 21 - 23/ ^^^Fe). Data: O - 
Engelmann et al. (1990), A - Binns et al. (1988). 
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Figure 3: Model interstellar and modulated (500 MV) ratios. Left: ^''Be/^Be. Data from Lukasiak et al. 
(1994a) (D - Voyager-1,2, o - IMP-7/8, A - ISEE-3) and Connell (1998) (• - Ulysses). Right: ^^Al/ ^^Al. 
Data: o - Lukasiak et al. (1994b), • - Simpson and Connell (1998). Note that the data points shown are at 
the measured (not interstellar) energies. 
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Figure 4: Model interstellar and modulated (500 MV) ratios. Left: ^''Cl/ CI. Data: Connell et al. (1998); 
Right: 5^Mn/ Mn. Data: • - Duvernois (1997), o - Lukasiak et al. (1995), A - Leske (1993). 



